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SUMMARY 


TWb  paper  suggests,  aa-.did  an  earlier  one,  [1]  that  points 
In  n-space  produced  by  congruential  random  number  generators  are 
too  regular  for  general  Monte  Carlo  use.  Regularity  was 
established  in-- fir]  for  multiplicative  congruential  generators  by 
showing  that  all  the  points  fall  in  sets  of  relatively  few  parallel 
hyperplanes.  The  existence  of  many  containing  sets  of  parallel 
hyperplanes  was  easily  established,  but  proof  that  the  number  of 
hyperplanes  was  small  required  a  result  of  Minkowski  from  the 
geometry  of  numbers — a  symmetric,  convex  set  of  volume  2n  must 
contain  at  least  two  points  with  integral  coordinates.  The  present 
paper  takes  a  different  approach  to  establishing  the  coarse  lattice 
structure  of  congruential  generators.  It  gives  a  simple,  self- 
contained  proof  that  points  in  n-space  produced  by  the  general 
congruential  generator  r^+^  =  +  b  mod  m  must  fall  on  a  lattice 

with  unit-cell  volume  at  least  m11  \  There  is  no  restriction  on  a 
or  b;  this  means  that  all  congruential  random  number  generators 
must  be  considered  unsatisfactory  in  terms  of  lattices  containing 
the  points  they  produce,  for  a  good  generator  of  random  integers 
should  have  an  n-lattice  with  unit-cell  volume  1. 


The  Lattice  of  a  Random  Number  Generator 


Suppose  we  define  the  n -lattice  of  a  random  number  generator 
as  follows:  if  the  generator  produces  integers  r^rj,^,...  let 

ti^  *  (ri,r2 » •  •  •  *rn) »  ^2  "  (r2  »r3 » •  •  •  »rn+i^  >  •  •  •  be  set  possible 

points  in  n-space  formed  from  n  successive  r's.  The  n-lattice  of 
the  generator  is  the  set  of  all  integral  linear  combinations  of 
points  from  this  set  translated  to  include  the  origin,  i.e.,  all 
integral  linear  combinations  of  the  points 

TT2'“TT^,Tr^“7T^,'IT^“7T^,»**  .  (1) 

The  unit-cell  volume  of  the  n-lattice  is  the  greatest  common 
livisor  of  the  volumes  of  parallelepipeds  formed  from  any  n+1 
points  of  the  lattice;  the  volume  of  such  a  parallelepiped  is  the 
determinant  with  rows  formed  by  subtracting  one  of  the  points  from 
each  of  the  other  n  points. 

The  unit-cell  volume  may  be  considered  a  generalization  of 
the  idea  of  the  greatest  common  divisor  of  a  set  of  zero-translated 
integers,  and  even  a  few  dozen  points  in  n-space  with  truly  random 
integer  coordinates  is  virtually  certain  to  have  an  n-lattice  with 
unit-cell  volume  1.  The  following  theorem  shows  that  the  lattice 
structure  of  every  congruential  random  number  generator  is  far  too 
gross  to  make  the  generator  suitable  for  general  Monte  Carlo  use: 
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THEOREM.  Let 

t1  -  (1»T  (1)  ,T2  (1) . Tn_1  (1)  ) 

t2  -  (2,T (2) ,T2 (2) . Tn"1(2)) 

Tm  -  (0,T(0),T2(0) . Tn_1  (0) ) 

be  the  set  of  all  possible  points  in  n-spaoe  whose  coordinates  are 
generated  successively  from  an  initial  coordinate  by  a  linear  trans¬ 
formation  T  on  the  ring  of  reduced  residues  of  some  modulus  m: 

T(x)  =  ax  +  b  mod  m  0  <_  T(x)  <  m, 

or3  using  the  greatest  integer  notation 3 

T(x)  -  ax  +  b  -  m[ (ax+b)/m] . 

Then  all  of  the  points  »T2» •  •  •  »Tm  on  a  lattice  with  unit-cell 
volume  mn  . 

The  proof  is  not  very  difficult  and  the  case  n  *  3  will  serve 
to  describe  the  general  situation.  Since  the  volume  of  the  unit  cell 
of  a  lattice  is  the  greatest  common  divisor  of  the  volumes  of 
parallelepipeds  formed  from  sets  of  n+1  points,  it  suffices  to 
prove  that,  for  any  reduced  residues  r,  s,  t,  v, 

r-v  T(r)-T (v)  T2(r)-T2(v) 

s-v  T(s)-T(v)  T2(s)-T2(v)  e  0  mod  m""1. 

t-v  T(t)-T(v)  T2  (t  )-T2  (v) 
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Now  it  is  easy  to  verify  that 

TJ(r)  -  Tj(v)  -  a^(r-v)  s  0  mod  m  (2) 

2 

(since,  for  example,  T^(r)  ■  a^r  +  ab  +  b  -  m[5J r+ab+b j  ^  ^  Subtracting 

m 

a'*'  times  the  first  column  from  the  i^  column,  for  i-2,3  will 

produce  a  determinant  whose  columns,  except  the  first,  have  m  as  a 

factor.  In  general,  then,  the  determinant  will  have  mn  ^  as  a  factor. 

Since  there  are  m  distinct  points  ,  the  unit-cell  volume  of 

1  m 

their  lattice  will  be  exactly  m11  1.  Points  produced  by  any  particular 
congruential  generator  will  be  a  subset  of  the  x's  and  will  have  a 
lattice  with  unit-cell  volume  at  least  mn  ^ . 

Note  that,  by  virtue  of  (2),  every  zero-translated  point  in  n-space 
will  have  the  form 

(x,ax-ym,a^x-zm, . . .) 

and  this  form  readily  provides  a  basis  for  the  lattice — for  example,  the 
rows  of  this  matrix  are  a  basis  of  the  4-lattice: 

f  1  a  a2 

0  m  0 

0  0m 
\0  0  0 


-4- 


Relation  to  Sets  of  Parallel  Hyperplanes. 

The  paper  cited  above  [1]  suggested  that  the  crystalline  structure 
of  multiplicative  congruential  generators  was  too  crude  by  showing  that 
points  in  n-space  produced  by  such  generators  must  lie  in  a  set  of  less 
than  (n!m)^n  parallel  hyperplanes.  The  above  theorem  shows  that  every 
congruential  generator  produces  points  in  n-space  which  fall  on  a  lattice 
with  unit-cell  volume  m11  To  relate  the  two,  imagine  the  "best  possible" 
lattice,  with  cubic  structure,  and  with  one  of  the  sets  of  parallel  faces 
perpendicular  to  the  longest  line  through  the  cube  of  points  with  integer 
coordinates  in  the  range  0  to  m.  The  length  of  the  diagonal  is  m/n, 
and  the  length  of  a  tide  of  the  cubic  unit-cell  is  (mn  Dividing 

the  length  of  the  diagonal,  mt/n,  by  the  distance  between  parallel 
hyperplanes,  m^n  1^n,  we  get  this  bound  for  the  number  of  hyperplanes 
containiug  all  the  points  of  a  congruential  random  number  generator: 

m1/,n.  This,  if  true,  would  be  an  improvement  on  the  previous 
bound,  (n!m)^ri.  Can  the  argument  be  made  rigorous?  The  question 
is  mainly  academic,  for  in  either  case  the  bound  is  too  low  to  make 
congruential  generators  suitable  for  general  use. 
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